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Abstract 

Communication systems in practice are subject 
to many technical/technological constraints and 
restrictions. MIMO processing in current wire¬ 
less communications, as an example, mostly em¬ 
ploys codebook based pre-coding to save compu¬ 
tational complexity at the transmitters and re¬ 
ceivers. In such cases, closed form expressions 
for capacity or bit-error probability are often 
unattainable; effects of realistic signal processing 
algorithms on the performance of practical com¬ 
munication systems rather have to be studied 
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in simulation environments. The Vienna LTE- 
A Uplink Simulator is a 3GPP LTE-A standard 
compliant link level simulator that is publicly 
available under an academic use license, facili¬ 
tating reproducible evaluations of signal process¬ 
ing algorithms and transceiver designs in wireless 
communications. This paper reviews research re¬ 
sults that have been obtained by means of the 
Vienna LTE-A Uplink Simulator, highlights the 
effects of Single Carrier Erequency Division Mul¬ 
tiplexing (as the distinguishing feature to LTE-A 
downlink), extends known link adaptation con¬ 
cepts to uplink transmission, shows the impli¬ 
cations of the uplink pilot pattern for gather¬ 
ing Channel State Information at the receiver 
and completes with possible future research di¬ 
rections. 

1 Introduction 

Current cellular wireless communications em¬ 
ploys Universal Mobile Telecommunications Sys¬ 
tem (UMTS) Long Term Evolution (LTE) as 
the high data rate standard [1]. The increas¬ 
ing demand of high data traffic in up- and 
downlink forces engineers to push the limits 
of LTE [2], e.g., through enhanced multi-user 
Multiple Input Multiple Output (MIMO) sup- 
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port m SI, Coordinated Multipoint (CoMP) 
transmission/reception [5l[6] as well as improved 
Channel State Information (CSI) feedback al¬ 
gorithms [7]. The authors of [8] predict fur¬ 
ther evolution of existing LTE/ LTE-Advanced 
(LTE-A) systems in parallel to development of 
new radio-access technologies operating at mil¬ 
limetre wave frequencies until 2020 and beyond. 
Fair comparison of novel signal processing al¬ 
gorithms and transceiver designs has to assure 
equal testing and evaluation conditions to enable 
reproducibility of results by independent groups 
of researchers and engineers [9]. For perform¬ 
ing system-level simulations m, [ii] or [12] are 
freely accessible options. For link level there are 
mainly commercial products available that facil¬ 
itate reproducible research, such as, is-wireless 
LTE PHY LAB [T3| or Mathworks LTE System 
Toolbox [nj. To the best of the authors knowl¬ 
edge, however, the Vienna LTE simulators are 
the only suite of simulation tools for LTE sys¬ 
tem and link level publicly available under an 
academic use licence, thus, free of charge for aca¬ 
demic researchers all over the world. In this pa¬ 
per, we introduce the latest member of the fam¬ 
ily of Vienna LTE simulators, that is, the Vienna 
LTE-A uplink link level simulator, and highlight 
our research conducted by means of this simula¬ 
tor. 

The outline of this article is as follows: We 
start with a brief re-capitulation of the LTE- 
A specifics and introduce the modulation and 
multiple access scheme and the employed MIMO 
signal processing of LTE-A uplink in Section [2l 
We then develop a matrix model describing the 
input-output relationship of the LTE-A uplink 
and present Signal to Interference and Noise Ra¬ 
tio (SINK) expressions for Single Carrier Fre¬ 
quency Division Multiplexing (SC-FDM) as well 
as Orthogonal Frequency Division Multiplexing 


(OFDM). The OFDM SINK expression and the 
performance of OFDM will serve as reference to 
study to effects of DFT-spreading imposed by 
SC-FDM. 

In Section [S] we investigate the physical layer 
performance of SC-FDM and OFDM, compar¬ 
ing Bit Error Ratio (BER) and Peak to Aver¬ 
age Power Ratio (PAPR). BER for LTE SISO 
transmissions were already analysed in link-level 
simulations by [laiiiiiT] and semi-analytically 
by [181119j . By means of our simulator, we re¬ 
produce these results and provide bounds to pre¬ 
dict the performance of SC-FDM with respect to 
OFDM. The insights gathered by the BER sim¬ 
ulations allow us to interpret the difference in 
throughput obtained by OFDM and SC-FDM, 
as discussed in Section 01 

Based on the SINR expressions developed in 
Section [2l we present a limited feedback strategy 
for link adaptation in Section 0] and contrast the 
performance of LTE uplink with channel capac¬ 
ity and other performance upper bounds that ac¬ 
count for practical design restrictions [20| . Until 
Section [5] we assume perfect CSI at the receiver. 
The remaining sections will describe methods to 
obtain CSI at the receiver. 

In Section [5l we highlight and describe the De¬ 
modulation Reference Signal (DMRS) structure 
employed in LTE-A uplink to facilitate channel 
estimation of the time-frequency selective wire¬ 
less channel. 

Based on the obtained insights, we elaborate 
on the basic concept of Discrete Fourier Trans¬ 
form (DFT)-based time domain channel estima¬ 
tion in Section [6] and review alternative code / 
frequency domain methods that can outperform 
DFT-based schemes m- 

Due to the increasing number of mobile users 
that stay connected while travelling in cars or 
(high speed) trains, we then shift our focus to 
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Resource Block 


high velocity scenarios. Such scenarios entail 
high temporal selectivity of the wireless channel, 
rendering accurate channel interpolation very 
important to sustain reasonable quality of ser¬ 
vice. We introduce and investigate basic con¬ 
cepts of channel interpolation in Section [71 

We briefly discuss open questions for future 
research in Section [8] and conclude in Section [9l 
Details to the handling of the simulator are pro¬ 
vided in p2] . 


Notation 

Matrices are denoted by bold upper-case letters 
such as H and vectors by bold lower-case let¬ 
ters such as h. The entries of vectors and ma¬ 
trices are accessed by brackets and subscripts, 
e.g., [h]k and [H]k^n- Spatial layers or receive 
antennas are denoted by superscripts in braces, 
e.g., The superscripts {■)'^ and (•)^ ex¬ 

press transposition and conjugate transposition. 
II • lb) II • Iloo and II • \\f symbolizes the Euclidean-, 
the Maximum- and the Frobeniusnorm, respec¬ 
tively. The entrywise (Hadamard) product is de¬ 
noted by 0 and the Kronecker product by 0. 
The all ones vector/matrix is denoted by 1. The 
operator X = Diag(a;) places the vector x on the 
main diagonal of X and conversely the operator 
X = diag(X) returns the vector x from the main 
diagonal of X. A block-wise Toeplitz (circulant, 
diagonal) matrix is a block matrix with each ma¬ 
trix of Toeplitz (circulant, diagonal) shape. The 
size of matrices is expressed via their subscripts, 
whenever necessary. 



Figure 1: The LTE-A uplink resource grid. 


2 LTE-specific System Model 
and SINK 

LTE operates on a time-frequency grid as shown 
in Figured! The number of subcarriers is always 
a multiple of twelve; twelve adjacent subcarriers 
over seven (or six ~ in case of extended Cyclic 
Prefix (CP) successive OFDM symbols are called 
Resource Block (RB). Each RB thus consists of 
12 X 7 (12 X 6) Resource Elements (REs), corre¬ 
sponding to the different time-frequency bins. A 
detailed description of LTE up- and downlink is 
available, e.g., in [23] . 

We focus one those details, necessary to de¬ 
scribe our system model at time rii|. LTE em¬ 
ploys OFDM (A) as physical layer modulation 
and multiple access scheme in the downlink and 
SC-FDM(A), i.e., DFT-spreaded OFDM, in the 
uplink. In a SC-FDM model, OFDM can be con¬ 
sidered a special case. The major difference is 
an additional spreading and de-spreading stage 

^Note that we use the symbol n as time index and the 
vector n for noise, the distinction should be clear from 
the context. 
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Figure 2: The LTE-A uplink transceiver. 


at the transmitter and receiver, highlighted via 
dashed boxes in Figure [2j The common parts of 
the system model will be described from left to 
right. 

Right after the DFT spreading, the DMRS 
are inserted. DMRS will be considered later for 
the purpose of Channel Estimation (CE). Next, 
MIMO precoding is carried out, exploiting a set 
of semi-unitary precoding matrices W, pooled 
in the precoder codebook W, as defined in [T]. 
For LTE-A uplink transmission, the precoding 
matrix applied for a given user is equal for all 
RB assigned to this user. In case of spatial mul¬ 
tiplexing, each spatial layer is transmitted with 
equal power. 

Each antenna is equipped with its own OEDM 
modulator, consisting of subcarrier mapping. In¬ 
verse East Fourier Transform (IFFT) and an 
CP addition. To cope with the channel dis¬ 
persion and to avoid Inter Symbol Interference 
(ISI), LTE employs a CP. As a result of mul¬ 


tipath propagation a previous symbol may over¬ 
lap with the present symbol, introducing ISI and 
impairing the orthogonality between subcarriers, 
i.e., causing Inter Carrier Interference (ICI) |24j . 
Normal and extended CP lengths, with a respec¬ 
tive duration of 4.7//S and 16.7//S, are standard¬ 
ized, enabling a simple trade-off between ISI im¬ 
munity and CP overhead. 

At the transmitter, processing occurs in 
reversed order. First the OFDM demodula¬ 
tion / FFT takes place to get back into the 
frequency domain. The immunity to multipath 
propagation (stemming from the CP) allows to 
employ one-tap frequency domain equalizers F 
without performance loss. At last, de-spreading 
delivers the data estimates. 

All this previously informally described process¬ 
ing is linear and we are able to formulate a 
matrix-vector input-output relationship between 
a (stacked) data-vector x and its estimate x. For 
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X = {II ® ^Nsc) ^ ® DNppj,PremCp) H (I^Vt ® PuddCP^NFFT^) ^Nsc) L ® X 

IlNsc ^lNsc + [II® Dnsc) F {INr® D^pp^Pj-eraCv) '<X 

for OFDM for OFDM '-v/-' 

/ ^ \ > -'''^-V n 

= {II ® Dnsc) FH^q {II ® Dnsc) x + fi 

=Kx + n= IQKx + {K — IQK)x + h . (1) 

desired signal infra- and interlayer interference 


simplicity we assume that the channel stays con¬ 
stant during one OFDM symbol. A detailed 
system description based on [25] can be found in 

m- 

In order to adapt the data transmission to 
the current channel state, LTE-A applies lim¬ 
ited feedback; a comprehensive specification fol¬ 
lows in Section 01 Limited feedback is depicted 
via the feedback arrow in Figure [2j The data 
vector G gf layer I G {1,..., L} con¬ 

tains modulated symbols for each of the Nsc 
subcarriers. The number of transmit layers de¬ 
pends on the LTE-A specific Rank Indicator (RI) 
feedback. The data symbols are coded with a 
punctured turbo code whose rate is determined 
by the Channel Quality Indicator (CQI). Sub¬ 
sequently, the codewords are mapped onto a 
Quadrature Amplitude Modulation (QAM) al¬ 
phabet (4/16/64 QAM), where the size of the al¬ 
phabet depends on the CQI as well. All x^’‘'> are 
stacked into one vector x G which 

layer-wise spreading and joint precoding - ac¬ 
cording to the Precoding Matrix Indicator (PMI) 
- of all subcarriers takes place. The subsequent 
OFDM modulator consists of the localized sub¬ 
carrier mapping M, mapping Nsc subcarriers 
to the center of an NppT point IFFT, and the 
addition of the CP. 


Depending on the level of abstraction, our sys¬ 
tem model can be described via different channel 
matrices. The physical baseband time domain 
channel is described by a block-wise Toplitz ma¬ 
trix H G , with 

Np transmit and Np receive antennas, which 
turns block-wise circulant (iTcir) after addition 
(-Paddcp) and removal (Premcp) of an appropri¬ 
ately chosen CP of length Ncp- Finally, it turns 
diagonal after the IFFT and FFT on the trans¬ 
mitter and receiver, respectively. An example of 
the Toplitz and diagonal structured channel is 
demonstrated in Figure [3] (a) and (b), respec¬ 
tively. 


P-diag = {lNii®DNpp^PremCp) H [lN^'S>PaddCpF>^^^^) 

( 2 ) 

The last step of the OFDM de-modulator is 
the reversal of the localized subcarrier mapping 
. The effective MIMO channel from L trans¬ 
mit layers to Np receive antennas, incorporating 
the precoder, the OFDM modulator, the time- 
domain MIMO channel H and the OFDM de¬ 
modulator, is abstracted to one block matrix 
iTeflf- This greatly facilitates the readability of 
all formulas later on. 
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(a) Magnitude values of the Toplitz 
structured MIMO time domain channel 
matrix with 2 transmit antennas and 4 
receive antennas. 
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(b) Diagonalized MIMO channel ma¬ 
trix after CP addition/removal and 
IFFT/FFT at the transmitter and re¬ 
ceiver, respectively. 
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(c) Input-Output matrix K employing 
an MMSE equalizer. The non-diagonal 
elements represent the intra- and inter¬ 
layer interference. 



( 3 ) 


The additive noise is assumed independent 
across antennas and is distributed zero mean, 
white Gaussian ~ CV{0,cr^/}, i G 

The stacked noise vector n = 
^ is thus zero mean, 

white Gaussian as well. 

The frequency domain one-tap equalizer F is 
chosen conforming to different criteria, either the 
Zero Forcing (ZF) criterion, which removes all 
channel distortions at risk of noise enhancement, 
or the Minimum Mean Squared Error (MMSE) 
criterion, that tries to minimize the effects of 
noise enhancement and channel distortion. 

After the de-spreading operation the data es¬ 
timates X of the noisy, received signal are given 
in Equation ([T]), with the before mentioned con¬ 
venient abbreviation ([3]) and DNppt is the DET 
matrix of size Nfft- ■ 


2.1 SC-FDM SINK 

The special structure of Equation ([T]), due to the 
frequency domain one tap equalizer and the DET 
spreading, yields a block-wise circulant input- 
output matrix, cf. Eigure [3] (c), 

K={Il® ® Dnsc) • ( 4 ) 

This block-wise circulant structure produces a 
constant post equalization and post spreading 
SINK over all subcarriers within one layer [26]. 
The detailed derivation is provided in the ap¬ 
pendix. 
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Eigure 3: Examples of different channel abstrac¬ 
tions. 
























with the selection vector 


SinrSC-fdm. (0 ^ ^5) 

ai||S«Ff/,ff|||-^|lT^^S«diag(FJ/eff)| +a2||S«F|||’ 

where 

5« = (0 Insc 0) , (6) 

selects that part of FH^s effecting the layer. 
The second moment of the zero mean symbols 
equals the baseband transmit power as LTE- 
A has standardized semi-unitary precoders W, 
so that the overall transmitter (spreading, pre¬ 
coding and OFDM modulation) is unitary. 

2.2 OFDM SINK 

In contrast to SC-FDM, no spreading takes place 
for OFDM. The dashed boxes in Figure [2] are 
replaced by identity matrices; they are simply 
omitted. Thus, different sub carriers k are or¬ 
thogonal/independent and the equalizer treats 
the corresponding subcarrier channel only. 
We use the subscript k to denote the relevant 
part of the full channel matrix iTeff for the 
subcarrier. The corresponding indices within the 
diagonal matrix iTdiag are InrxNt ® Diag (e^), 
with the canonical base vectors e^. Using this 
notation, the effective subcarrier channel iTfc G 

C^rxL ig 


s^^'> = ( 0...0 1 0...0) , (9) 

with appropriate number of zeros and a one at 
the position. 


3 SC-FDM Effects 


3.1 Peak to Average Povi^er Ratio 

SC-FDM is employed as the physical layer mod¬ 
ulation scheme for LTE uplink transmission, due 
to its lower PAPR compared to OFDM [27]. 
Lower PAPR, or similarly lower crest factor, 
leads to reduced linearity requirements for the 
power amplifiers and to relaxed resolution spec¬ 
ifications for the digital-to-analog converters at 
the user equipments, entailing higher power effi¬ 
ciency. 

The Vienna LTE-A uplink simulator calcu¬ 
lates the discrete-time baseband PAPR with the 
default oversampling factor J = 4 [28] . The 
discrete time signal on transmit antenna t G 
{!,..., Nt} is therefore calculated as 




1 




Njpp'j' — 1 


(t) 

pre 


• 27rmfc 

o'' JNj^r^rp 


- k 


fft 

0 < m < — 1 , 


( 10 ) 


and Ffc is its linear one tap equalizer. The SINR 
formula is quite similar to the SC-FDM case, ex¬ 
cept that the SINR shows subcarrier dependency 
now. The SINR vector at layer I reads 


0-2 |s(')diag(FfcJi'fc)| 


- k 
2 


( 8 ) 


tT2||s(0Ffc//fc|li-fT2|smdiag(Ffc//D| +fT2||s(0Ffc| 


where is the transmit vector right after 

precoding and before the IFFT at transmit an¬ 
tenna t. The PAPR of the stacked vector Stx = 


((■ 




tx 


>Kx 


(Nrp) 


)^) is calculated as 


PAPR{stJ 


l<t<JVj. 0<m<JNppj,-l 


„(‘)1 




(11) 


NTNpp^Wdiag (stxsS) || Oo/||^tx||2 5 
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Figure 4: PAPR for SC-FDM and OFDM for 
different bandwidths (1.4 MHz and 10 MHz) and 
modulation alphabets (4/64 QAM). 


where the Euclidean norm in the denominator 
serves as an estimate for the ensemble average. 

Figure 0] depicts the PAPR of OFDM and 
SC-FDM obtained for different system band- 
widths. Already for a small bandwidth (1.4 
MHz), there is a significant reduction for SC- 
FDM over OFDM. With increasing bandwidth 
OFDM’s PAPR grows and the gains obtained 
by SC-FDM become more and more pronounced. 
The PAPR also depends on the modulation al¬ 
phabet; the smaller the alphabet, the smaller the 
PAPR. This effect is illustrated in dotted lines in 
Figured! where we have shown the PAPR of 4- 
QAM, exemplarily. 

3.2 BER Comparison over Frequency 
Selective Channels 

The additional spreading of SC-FDM leads to an 
SINR expression that is constant on all subcar¬ 
riers as for single carrier transmission, legitimat¬ 
ing its name. The aim of this subsection is to 
analyse the SINR expression more carefully for 


the Single Input Single Output (SISO) cas^ and 
draw conclusions on BER performance. 

We focus on the two most prominent equal¬ 
izer concepts and start with the ZE equalizer, for 
whom the SC-EDM Signal to Noise Ratio (SNR) 
expression ([5|) reduces to the harmonic mean 


SNR 


SC-FDM 

ZF 




1 


^ Nsc 

^SC \Hk\ 


( 12 ) 


whereas the OEDM expression ([8|) is sub-carrier 
dependent and becomes proportional to the 
channel transfer function 

'^n 

The average OEDM SNR 


snr9|°^ 


O'x 1 

alNsc 


^sc 


k=l 


(14) 


yields an upper bound on the Single Carrier Ere- 
quency Division Multiple Access (SC-EDMA) 
SNR due to the harmonic mean - arithmetic 
mean inequality m- 

SNR|^“™^ < SNRgI™ (15) 


Equality in Equation (I15p holds if and only if the 
channel is frequency flat. The difference between 
the harmonic mean and the arithmetic mean 
gets increasingly pronounced, the more selective 
the channel becomes. We therefore expect the 
(uncoded) BER of SC-EDM and ZE equaliza¬ 
tion to perform worse than OEDM, which is 
also validated by simulations. The BER sim¬ 
ulations were carried out with CQI = 4 on a 

^The reduction to SISO is done to make our results 
comparable even to older Frequency Domain Equalization 
(FDE) works, e.g., [29] 
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PedB channel m- This Modulation and Coding 
Scheme (MCS) employs 4-QAM and has an effec¬ 
tive code-rate of 0.3008. As expected, the BER 
performance of SC-FDM is worse than OFDM, 
both shown in Figure [S] (a) in solid lines. Due 
to the spreading SC-FDM already expends all 
channel diversity and coding does not increase 
the SNR slope of the BER curve. This manifests 
in an almost parallel shift of the BER curve for 
SC-FDM, as visual in Figure [5] (a) in red dashed 
lines. None exploited diversity allows OFDM to 
increase the BER slope considerably, cf. Figure 
[5] (a) blue dashed line. 

The MMSE SINR expression is less intuitive 
and for the purpose of comparison, similar math¬ 
ematical transformations as in |32] and |19j are 
required to arrive at 



(b) MMSE receiver 


Figure 5: BER comparison between OFDM and 
SC-FDM for a SISO PedB channel with 5 MHz 
bandwidth and fixed CQI = 4 transmission. 


SINR 


SC-FDM 

MMSE 



1-; 


.2 , ^SC 

n 1 

?T Nor ^ 

^ fe=l 



1 

^SC 


^SC 


E 


k=l 



(16) 



1 


, ^sc 




The detailed derivation is shown in the appendix. 
The denominator of Equation (11611 is regularized 
and less sensitive to spectral notches. 

An upper bound on the SINR can be obtained 
via the maximum of the transfer function Hj, 


SINR 


SC-FDM 

MMSE 



= ^maxfc|iTfcp. (17) 

^ n 


In the low SNR regime ^ this bound 

2 

becomes tight. The higher the inverse SNR % 
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in relation to the maximum of the transfer func¬ 
tion, the tighter the bound becomes. As Equa¬ 
tion (I14p - the average OFDM SNR - can never 
be larger to its maximum entry (only equal for 
frequency flat channels), the SC-FDM SINK ex¬ 
pression (fT6]l is larger than at low SNR, 

as predicted by the tight bound (fTTl) : a lower 
BER is expected. Again, this presumption is val¬ 
idated by our simulation, showing that the un¬ 
coded BER is lower for SC-FDM as for ODFM, 
cf. Figure [5] (b) in solid lines. Although the 
uncoded BER shows superior performance, the 
coded BER is lower for OFDM due to the cod¬ 
ing gains stemming from channel diversity, cf.. 
Figure [5] (b) dashed lines. 

A bound for the Maximum Likelihood (ML) 
detection performance was derived in [33]. As 
bandwidth increases the slope of the BER curve 
achieved with MMSE receivers tends to the slope 
of ML detection, demonstrating the full exploita¬ 
tion of channel diversity by the MMSE equalizer, 
cf., Figure [S](b) black line. 

4 Link Adaptation 

In this section, we first investigate the through¬ 
put performance of LTE-A uplink employing 
ideal rate adaptation and compare SC-FDM 
transmission to OFDM with ZF and MMSE 
receivers. Then, we extend our single-user 
MIMO CSI feedback algorithms proposed for 
LTE downlink in [3Tj to LTE uplink and evaluate 
their performance comparing to the throughput 
bounds developed in [20| . We also highlight some 
important basic differences between link adapta¬ 
tion in LTE up- and downlink transmissions. 

®For the SISO case, ZF and MMSE equalizers perform 
equivalent for OFDM and 4-QAM. 


4.1 Performance with Ideal Rate 
Adaptation 


As demonstrated in the previous section, SC- 
FDM provides a signihcant advantage in terms 
of PARR over OFDM, thus relaxing linear¬ 
ity requirements of radio frequency power am¬ 
plifiers for user equipments. Yet, this comes 
at the cost of coded BER degradation since 
channel diversity is lost and the performance 
is mostly dominated by the weakest subcar¬ 
rier of a user, especially with ZF receivers; 
c.L, (I12p . This diversity loss cannot be recov¬ 
ered by forward-error-correction channel coding, 
since the DFT-spreading applied with SC-FDM 
effectively causes an averaging over SINR ob¬ 
served on all scheduled subcarriers according 
to (HD. As a consequence, SC-FDM transmis¬ 
sion over frequency selective channels achieves 
worse throughput than OFDM. This is demon¬ 
strated in Figure (G] where we cross-compare 
the achievable rate, as defined in Equation (fTSP 
and (USD, and the actual throughput of SC- 
FDM and OFDM transmission as obtained by 
the Vienna LTE-A uplink simulator. We con¬ 
sider single-user transmission over 5 MHz band¬ 
width assuming = Nji = 2 antennas at 
the user and the base station and L = 2 spa¬ 
tial layers. The precoder is selected as a scaled 
identity matrix: W = I/^/LIl. We con¬ 
sider transmission over independent and iden¬ 
tically distributed frequency-selective Rayleigh 
fading channels, emphasizing the difference be¬ 
tween OFDM and SC-FDM. The achievable rate 
in bits per OFDM/SC-FDM symbol with Gaus¬ 
sian signalling and equal power allocation over 
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(a) ZF receivers 



Figure 6: Throughput comparison of OFDM and SC-FDM with rate adaptation and 2x2 Rayleigh 
fading channels of 5 MHz bandwidth. 


subcarriers and spatial layers is calculated as 
Nsc L 

rOFDM = £ ^ iog2 (l + [sINRO™^’ « 


k=l1=1 


- 

(18) 


R 


SC-FDM 


l=l 


Nsc logs (l + SINRSC-™^’ 1^1 

(19) 


with the receiver-specific post-de-spreading 
(post-equalization) SINRs from ([2]) and ([5]), re¬ 
spectively. 

We observe a significant loss of achievable rate 
of SC-FDM transmission compared to OFDM in 
Figure m which is especially pronounced with ZF 
receivers due to noise enhancement. In Figure [H 
we also show the actual rate achieved by LTE 
uplink SC-FDM transmission with ideal rate 
adaptation and compare to the performance ob¬ 
tained by OFDM transmission; the correspond¬ 
ing curves are denoted by LTE rate. We deter¬ 
mine the performance of ideal rate adaptation by 


simulating all possible transmission rates, corre¬ 
sponding to CQIl to CQI15, and selecting at 
each subframe the largest rate that achieves er¬ 
ror free transmission. The figure also shows the 
throughput of the individual CQIs. We observe 
a gap between the LTE throughput with OFDM 
and SC-FDM that is similar to the gap in terms 
of achievable rate. Notice that the performance 
loss with MMSE receivers is significantly smaller 
than with ZF detection, since MMSE avoids ex¬ 
cessive noise enhancement. 


We also observe in Figure 6(a) that the gain 
achieved by instantaneous rate adaptation, as 
compared to rate adaptation based on the long¬ 
term average SNR, is much larger for ZF SC- 
FDM than for ZF OFDM; this is evident from 
the distance between the curves with rate adap¬ 
tation {LTE rate) and the curves with fixed CQI. 
The reason for this behaviour is that the SNR 
of ZF SC-FDM shows strong variability around 
its means, since it is dominated by the worst- 
case per-subcarrier SNR according to (fT^ : the 
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average SNR over subcarriers of ZF OFDM, 
however, approximately coincides with its mean 
value. This implies that the optimal CQI of ZF 
SC-FDM can vary significantly in-between sub- 
frames, as reflected by the large average SNR 
variation required to increase the rate with fixed 
CQI from zero to its respective maximum. Yet, 
for ZF OFDM the throughput of the individual 
CQIs follows almost a step function; hence, rate 
adaptation can be based on the long-term aver¬ 
age SNR without substantial performance degra¬ 
dation!^ 

In case Nji > L, we can easily estimate the 
achievable rate of SC-FDM transmission: The 
per-layer SNR with ZF receivers is governed by 
the harmonic mean of the channel responses on 
the individual subcarriers, similar to (11211 


Letting iVgc —>• oo, we can replace the term in 
the denominator of (f20]l with its expected value 

— (23) 

k=l 

This expected value only exists in case Nr > 
L [36]. For Nr = L, the diagonal elements of H 
follow a heavy-tailed inverted Gamma distribu¬ 
tion [371 [38] with non-finite first moment. Yet, 
for Nr > L, which is a common situation in cel¬ 
lular networks since the base station is mostly 
equipped with far more antennas than the users, 
the expected value is 


SNR 


SC-FDM, (/) 
ZF 


1 

Nsc N,k=i 


({HkWf {HkW)y^ 

^ ^ \ I,I 

( 20 ) 


with iTfc G denoting the OFDM chan¬ 

nel matrix on subcarrier k. Assuming constant 
precoding and semi-correlated Rayleigh fading 


Hk 


HkC^, 



-CAT {0,1}, 


( 21 ) 


with Ct ^ (^NtxNt determining the spatial cor¬ 
relation at the user equipment side, the matrix 
in the denominator of (|20p follows a complex in¬ 
verse Wishart distribution with Nr degrees of 
freedom and scale matrix C = (VY^CtVF) 


H = {{HkW f {HkW)) ' ~ CWl^ {Nr, C} . 

( 22 ) 

■^Notice, however, that instantaneous rate adaptation 
for ZF OFDM can be advantageous in case of frequency- 
correlated channels [35]. 


Hence, we can estimate the achievable rate 
of SC-FDMA transmission over semi-correlated 
Rayleigh fading channels 


R 


SC-FDM 


IVSC ^ log2 

l=l 




\C]u 


(25) 


NscL 

log2 

^ 2/2 

+ log2 {Nr — L) 



(nti ic,,,) ' ) 



(26) 


Here ()26p resembles the high SNR approximation 
of the achievable rate of OFDM transmission 
with ZF detection as proposed in [39l Eq. (14)]; 
even more, for fixed L and letting Nr grow to 
infinity, (1^ and |39l Eq. (14)] tend to the same 
limit, due to channel hardening on each subcar¬ 
rier with growing number of receive antennas. 

In Eigure [T] we investigate the performance 
of the rate estimate (|25l) for Nr = L = 4 and 
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(a) Comparison of achievable rate and the rate esti- (b) Relative achievable rate of SC-FDM with respect 
mate of Equation (1251) . to OFDM. 

Figure 7: Achievable rate of OFDM and SC-FDM with ZF equalizers and growing number of 
receive antennas at fixed number of streams L = 4. 


varying number of receive antennas. We assume 
W = 1//L II and 

■ 1 0.9 ... 0.9' 



: 0.9 

0.9 ... 0.9 1 

and consider the smallest LTE bandwidth of 
-^sc = 72 subcarriers. We observe that the pro¬ 
posed estimate performs very well even at this 
small bandwidth; notice, though, that a more re¬ 
alistic channel model with correlation over sub¬ 
carriers may require larger bandwidth to validate 
the proposed estimate. Figure [7] also confirms 
the observation that single-user MIMO OFDM 
and SC-FDM with ZF detectors tend to the same 
limiting performance with increasing number of 
receive antennas. 

This statement, however, will not hold true if 
the total number of layers grows proportionally 


with the number of receive antennas. For ex¬ 
ample, multi-user MIMO transmission with ZF 
equalization and single antenna users achieves 
only a diversity order of Nr — L-\-1 [ID] , with L 
denoting the total number of layers being equal 
to the number of spatially multiplexed users. 
Hence, if L scales proportionally with Nr, chan¬ 
nel hardening on each subcarrier will not occur 
and thus the performance of OFDM and SC- 
FDM will not coincide. 

4.2 Performance with Realistic Link 
Adaptation 

Instantaneous rate adaptation is an important 
tool for exploiting diversity of the wireless chan¬ 
nel in LTE, by adjusting the transmission rate 
according to the current channel quality experi¬ 
enced by a user. LTE specifies a set of fifteen 
different MCSs; the selected MCS is signalled by 
the CQI. 
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(a) Absolute throughput 



Figure 8: Absolute and relative throughput of LTE uplink transmission over 4x4 VehA channels 
of 1.4 MHz bandwidth employing rate adaptation. We compare the performance of fixed rank, rank 
adaptive and PMI + rank adaptive transmission to the performance bounds proposed in |20j . 


LTE additionally supports spatial link adapta¬ 
tion by means of codebook based precoding with 
variable transmission rank. With this method, 
the precoding matrix W G (£^NtxL satisfying 
W^W = 1/L is selected from a standard de¬ 
fined codebook Wl of scaled semi-unitary matri¬ 
ces; furthermore, the number of spatial layers L 
can be adjusted to achieve a favourable trade-off 
between beamforming and spatial multiplexing. 
The selected precoder and transmission rank are 
signalled, employing the PMI and the RI. In 
single-user MIMO LTE uplink transmission, the 
same precoder is applied on all RBs that are 
assigned to a specific user, whereas frequency- 
selective precoding is supported in LTE down¬ 
link. 

There is a basic difference between the utiliza¬ 
tion of CQI, PMI and RI in up- and downlink 
directions of Frequency Division Duplex (FDD) 
systems. In downlink, the base station is reliant 
on CSI feedback from the users for link adapta¬ 


tion and multi-user scheduling m, since chan¬ 
nel reciprocity cannot be exploited in FDD. CQI, 
PMI and RI can be employed to convey such CSI 
from the users to the base station via dedicated 
feedback channels [35]. In the uplink, on the 
other hand, the base station can by itself de¬ 
termine CSI exploiting the Sounding Reference 
Signals (SRSs) transmitted by the users. In this 
case, CQI, PMI and RI are employed by the base 
station to convey to the users its decision on link 
adaptation that has to be applied by the users 
during uplink transmission. 

In principal, link adaptation must be jointly 
optimized with multi-user scheduling to opti¬ 
mize the performance of the system, since the 
effective SC-FDM SINR (and thus the rate) of 
a user depends on the assigned RBs according 
to ([5|). For reasons of computational complexity, 
however, we assume that the multi-user sched¬ 
ule is already fixed and determine link adapta¬ 
tion parameters based on this resource alloca- 
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tion. We modify the approach proposed in [M] do not consider signalling delays between the 
for LTE downlink transmission to determine the base station and the user. We employ the VehA 
link adaptation parameters in four steps: channel model m and compare the absolute and 

relative (to channel capacity) throughput to the 
1. Determine the optimal precoder for each performance bounds proposed in Chan- 


transmission rank L < min [Nt-, Nji) by max¬ 
imizing transmission rate 

L 

W{L) = argmax V / (sINRSC-fd^, (i) 
WgWl ^ 

(27) 

Here, function /(•) maps SINK to rate; this 
could be either an analytical mapping, such 
as (USD, or a mapping table representing the 
actual performance of LTE. In our simula¬ 
tions, we employ the Bit-Interleaved Coded- 
Modulation (BICM) capacity as proposed 
in [21], since LTE is based on a BICM ar¬ 
chitecture. 

2. Determine the optimal LTE transmission 
rates per layer for each L and W{L). We 
employ a target Block Error Ratio (BLER) 
mapping in our simulations to determine the 
highest rate that achieves BLER <0.1. 

3. Select the transmission rank L that maxi¬ 
mizes the sum rate over spatial layers, uti¬ 
lizing the LTE transmission rates determined 
above. 

4. Set the RI and PMI according to L and 
W{L), respectively and set the CQIs con¬ 
forming to the corresponding LTE transmis¬ 
sion rates. 


nel capacity is obtained by applying Singular 
Value Decomposition (SVD)-based transceivers 
and water-filling power allocation over subcarri¬ 
ers and spatial streams. Notice that we do not 
account for guard band and CP overheads when 
calculating the channel capacity; that is, we only 
consider subcarriers that are available for data 
transmission. The achievable channel eapaeity 
takes overhead for pilot symbols (DMRS and 
SRS) into account, corresponding to a loss of 
16.7% in our simulation. The achievable BICM 
bound additionally accounts for equal power allo¬ 
cation, codebook-based precoding, ZF detection 
as well as the applied BICM architecture as de¬ 
tailed in pU] . 

The performance of LTE uplink transmission 
with full link adaptation [PMI and rank adap¬ 
tive) is similar to the aehievable BICM bound 
but shifted by approximately 3 dB. Notice that 
the saturation value is not the same because the 
highest CQI of LTE achieves 5.55 bit/channel 
use, whereas the BICM bound saturates at 
6 bit/channel use. We also show the performance 
of LTE uplink when restricted to fixed precod¬ 
ing [rank adaptive) and fixed rank transmission 
[rank 1, 2, 3, 4)- We observe that rank adaptive 
transmission even outperforms the envelope of 
the fixed rank transmission curves, since instan¬ 
taneous rank adaptation selects the optimal rank 
in each subframe independently. In terms of rel¬ 
ative throughput, we see that LTE uplink with 


In Figure (8] we evaluate the performance of 
single-user MIMO LTE uplink transmission over 
Nt = Npi = 4 antennas with link adaptation, 
1.4 MHz system bandwidth and ZF receiver. We 


“Notice that the simulation setup is the same as em¬ 
ployed in [20] for the investigation of LTE downlink trans¬ 
mission, thus, facilitating the comparison of up- and 
downlink performance. 
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Figure 9: The LTE-A uplink reference symbol 
allocation in two slots (one subframe). 


ZF receivers achieves around 40-50 % of channel 
capacity; remember, though, that this does not 
include CP and guard band overheads. 


5 Reference Symbols 

In LTE uplink two types of reference signals are 
standardized. Eor CE and coherent detection, 
DMRS are exploited, while SRS are employed 
for channel sounding to enable frequency selec¬ 
tive scheduling. Eor the purpose of CE we will 
consider DMRS only. The reference symbols are 
defined in [ 1 ] and are explained in more detail 
in [32l 03]. As shown in Figur® DMRS are 
multiplexed in the resource grid at OFDM sym¬ 
bol time n = 3 in every slot. In a Physical Up¬ 
link Shared Channel (PUSCH) transmission of 
the LTE-A uplink, a DMRS occupies all sched¬ 
uled subcarriers. We assume that the user is 
assigned all Nsc subcarriers starting at 0 , i.e., 
/c E {0, 1,..., Nsc — !}• We denote the Zadoff- 
Chu (ZC) base sequence on Nsc subcarriers for 
one slot by r E The base sequences r 


are complex exponential sequences lying on the 
unit circle fulfilling 


I[^^]a:I = 1 - (28) 

In LTE-A the DMRS of different transmission 
layers in the same slot are orthogonal in terms 
of Frequency Domain Code Division Multiplex¬ 
ing (FD-CDM) [l2|. This is obtained by cycli¬ 
cally shifting the base sequence. Similar to |44] . 
DMRS on layer I for one slot are given by 

= Diag (r(')) = T^'^Diag (f) , (29) 

with the cyclic shift operator 

= Diag ..., ..., ^ 

(30) 

and the layer dependent cyclic shift ai. We 
further conclude from (I28]l - p0]l that = 

which implies = Insc- 

ploiting (f28]l . the product of two DMRS from 
layers I and u with l,u £ {1,..., L}, becomes 

= (T('))^T(“)Diag (f)^Diag (r) 

= Diag ... ... j ^( 31 ) 

with Aa = au — Oil being the cyclic phase shift 
between DMRS of two different spatial layers. 

The FD-CDM orthogonality can therefore be ex¬ 
ploited as 

^ ^ foruT^; 

(32) 

After transmission over a frequency selective 
channel, this orthogonality has to be exploited 
to separate all effective MIMO channels at the 
receiver. 


16 















6 Channel Estimation 

For channel estimation we exploit the system 
model only at symbol times, where reference sig¬ 
nals are allocated. For normal CP length this is 
the 4**^ symbol in each slot, i.e., n = 3 as shown 
in Figured Since we estimate the channel only 
at this single symbol time per slot, interpolation 
in time has to be carried out to obtain channel 
estimates for the whole resource grid. The effects 
of interpolation will be studied in Section [71 As 
illustrated in Figure O the DMRS are added af¬ 
ter DFT spreading, right before precoding. As 
the channel estimation takes place after the re¬ 
ceiver’s DFT, just before equalization, the sys¬ 
tem model for CE amounts to an OFDM system. 
The system model ([T]) therefore reads as 

y = + n' , (33) 

with (pre-equalization) noise 

n' = {l Nr® Af^DTVFFT-PremCp) , (34) 

and the stacked vector r consisting of DMRS 
^{i) g (j^ATscxi from all active spatial layers I E 
{!,...,L}, i.e. r = ((r(i))'^,...,(r(^))^)^. To 
consider the received signal separately for each 
receive antenna i, we can select the according 
part from y by left multiplying with the selector 
matrix from ([6]). The received signal = 
S^'^'^y on antenna i is given by 

= («!»’.■ ■ ■.«!»''’) ’■ + (35) 

1=1 

with the pre-equalization noise on 

receive antenna i and 

being the (i,/)'’'^ block of iTefr- Since 


is diagonal, we exploit the relations = 

Diag (r(^)) and = diag (iT^g^^) to estimate 
a channel vector rather than a matrix and rear¬ 
range terms in (1351) leading to 

yW (3g) 

1=1 

'-V-' 

R 

(i) 

with the stacked vector = 

..., of all effective chan¬ 

nels from L active layers to receive antenna i for 
which we will drop the subscript in the following. 

6.1 Minimum Mean Square Error Es¬ 
timation 

First we present a MMSE estimator where we 
exploit (j36p and estimate the stacked vector 
consisting of effective channels from all L active 
layers to receive antenna i. The MMSE CE for 
receive antenna i is given by 

^MMSE = argminEjll^^*^ - , (37) 

which leads to the well-known solution |45j 

^Smse = , 

(38) 

with C^d) = 

6.2 Correlation Based Estimation 

As a low complexity approach, we correlate 
(matched filter) the received signal with the ref- 
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(a) Mean Squared Error (MSE) curves of the pro¬ 
posed CE methods for a 2 x 2 transmission with L = 2 
on a Typical Urban (TU) channel. 



(b) Resulting uncoded BER curves for the proposed 
CE schemes for 4QAM. 



0 5 10 15 20 25 30 

SNR 


‘I’ D5% confidence 


Figure 10: Channel estimation performance comparison for block fading. 


erence symbol of layer I to obtain a channel es¬ 
timate for the effective channel from layer 
I to receive antenna i 

^ y{i) . ( 39 ) 

Inserting our system model (fObj) and exploit¬ 
ing dni), we obtain 

= ( 

U=1 

U^l 

inter-layer interference 

Here has the same distribution as 
since (Rm)H 

is unitary and introduces phase 
changes only, cf. ([2^ . Due to the allocation 
of DMRS on the same time and frequency 

resources on different spatial layers, the initial 

- {i,l) 

estimate h ’ of one effective MIMO channel 
actually consists of a superposition of all L 


effective MIMO channels to receive antenna i. 
The unintentional contributions in (jdOl) . from 
layers u ^ I are inter-layer interference, making 
it unsuited as initial estimate for coherent 
detection. Different methods to separate the 
different effective MIMO channels in (1401) will 
be presented in the following. 


6.2.1 DFT based Channel Estimation 

A well known approach for CE in LTE-A uplink 
is DFT based estimation |43] , which aims to sep¬ 
arate the MIMO channels contributing to (1400 
in time domain. For this the individual cyclic 
shift of each DMRS is exploited. Applying a 
DFT on the receive signal, the individual phase 
shifts will translate into shifts in time domain. 
This makes a separation of Channel Impulse 
Responses (CIR)s from different MIMO chan¬ 
nels possible by windowing. In our simulator we 
implemented a DFT based estimator as in [46] 


j,{u) ^(i) 
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ing factor 7 , is given by 


6.2.2 Averaging 


= {InscL + AQ^g) ' {R^^Yv 


(b 




(42) 


For physically meaningful channels, neighbour¬ 
ing subcarriers will be correlated within the co¬ 
herence bandwidth m- We utilize this prop¬ 
erty and exploit the DMRS structure to per¬ 
form frequency domain CE. As explained in EH, 
applying a sliding averaging on the initial esti- 
mate h from (j40p over 7 adjacent sub carriers 
(7 equals 1,2,4,4 for L equals 1,2,3,4, respec¬ 
tively) cancels the inter-layer interference, as¬ 
suming the channel to be frequency flat on these 
7 consecutive subcarriers. The sliding average is 
given by 


'^SAV 


t+7-1 


- k 




E E [■ 

t=/c—7+1 j=t 


h 




j j 


(41) 


for 7 < A: < Asc — 7 -|- 1. The second sum de¬ 
scribes the averaging of 7 elements while the first 
sum describes the shift of this averaging window. 


6.2.3 Quadratic Smoothing 

Another method exploiting channel correlations 
to estimate the channel in frequency domain is 
Quadratic Smoothing (QS). This scheme can¬ 
not remove the inter-layer interference entirely, 
which manifests in a higher error floor, but shows 
improved performance at lower SNR in return. 
As explained in m this estimation method, ex¬ 
ploiting the smoothing matrix Q and a smooth- 


Similar to (j4ip this can be interpreted as another 
way to cope with the inter-layer interference 
in (linp by post processing. This method does 
not use the DMRS structure explicitly but 
suppresses the interference by smoothing. It 
is therefore not able to cancel the complete 
inter-layer interference but shows a improved 
performance at low SNR. 


6.3 MSE and BER comparison 


We assume a single user 2x2 MIMO transmis¬ 
sion with Age = 72 subarriers, a fixed number 
of layers L = 2 and a TU channel model m 
at zero speed. We perform a simulation with 
one point extrapolation, cf.Section [3 and show 
the MSE curves of the proposed estimators in 
Figure 10(a) The DFT based CE (D-bCE) 


shows the highest error flow of all estimators 
at high SNR while the MMSE estimator, 
of course shows best performance over the 
whole SNR range. Compared to these two 
methods, the Sliding-Averaging estimator (jlTp . 
denoted by SAV, encounters a 8 dB SNR penalty 
when compared to MMSE, but comes closest 
to MMSE performance at high SNR. The 
quadratic smoothing estimation is denoted by 
QS and shows a significant improvement for low 
SNRs because it smooths over several observed 
channel coefficients. Quadratic smoothing 
performs uniformly better than D-bCE over the 
whole SNR range and comes close to 4dB to 
MMSE at low SNR. The high error floor shows 
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that QS is not able to cancel all the inter-layer 
interference. 

In terms of BER performance, at high SNR, nat¬ 
urally the estimation method with lowest MSE 
leads to the smallest BER. At low SNR, the 
difference in CE MSE translates into very small 
differences in BER, meaning we cannot gain too 
much from a good low SNR MSE performance 
of QS or MMSE estimation. Considering 
estimation complexity and that MMSE as well 
as QS require prior channel knowledge, SAV 
estimation is a good complexity performance 
trade-off. 


7 Channel Interpolation 

Under fast fading conditions additional effects 
influence the performance of LTE uplink trans¬ 
missions. Doppler shifts degrade the SINR by 
introducing velocity dependent ICI [38] whereas 
the SINR increases with increasing subcarrier 
spacing. The subcarrier spacing of 15 kHz 
that is used in LTE makes transmissions quite 
robust against ICI. The impact of ICI becomes 
only evident at high velocities and high SNR. 
Fig. mi (b) shows the BER for the case of perfect 
channel knowledge where the performance is 
only degraded by noise and ICI. At 200 km/h 
the BER saturates due to ICI at high SNR 
whereas ICI mitigation techniques [IS] show 
promising results to reduce this impact of ICI. 

Another effect that hampers LTE transmis¬ 
sions at high velocities are temporal channel in¬ 
terpolation errors. While in the LTE down¬ 
link the pattern used to multiplex data and 
reference symbols is a good trad-off between a 
small temporal and spectral spacing account- 


subframe u-1 subframe u , , , , 

□ data symbol 

■ pilot symbol 


time 1. slot 2. slot 

(a) Resource grid 




(b) 1 point extrapolation 



(c) 2 point linear interpolation 



(d) 3 point linear interpolation 



(e) 3 point spline interpolation 

Eigure 11: Channel interpolation techniques us¬ 
ing the estimates from (b) the actual slot, (c) 
the actual subframe and (d)-(e) the actual and 
previous subframe. 


ing for highly frequency selective channels and 
fast fading channels and a rather small over¬ 
head, this is different in the uplink. As shown 
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in Fig. [TT] (a) uplink DMRSs occupy the whole 
subband. While there is no need for interpo¬ 
lation over frequency, the temporal spacing is 
about twice the spacing of the reference sym¬ 
bols in the downlink. Furthermore, if frequency 
hopping is performed the number of adjacent pi¬ 
lots transmitted in the same subband is two for 
inter-subframe frequency hopping and only one 
for intra-subframe frequency hopping where fre¬ 
quency hopping is performed on a per-slot basis. 
Due to this special structure channel interpola¬ 
tion in the LTE uplink is a challenging prob¬ 
lem. Therefore we investigated various channel 
interpolation techniques using a single, two or 
three consecutive pilot symbols. Fig. [TT] (b)-(e) 
illustrates the channel interpolation techniques 
considered. The highest channel interpolation 
errors (Fig. (T^j (a)) are observed for 1 point ex¬ 
trapolation where the channel estimate obtained 
in a certain slot is used to equalize the sym¬ 
bols within that slot and no interpolation is per¬ 
formed at all. The higher the number of pi¬ 
lots involved in channel interpolation the lower 
the MSE gets. The results in terms of BER in 
Fig. [I2] (b) show a similar behaviour. 

For a measurement based comparison of inter¬ 
polation techniques using channel estimates form 
both, the previous and the subsequent subframe 
the reader is referred to [50] . 

8 Future Research Questions 

Until now our research efforts on the Vienna 
LTE-A Uplink Simulator have been concentrated 
on single links between user and base station, 
focusing on basic transceiver issues such as link 
adaptation and channel estimation. In future, 
our scope will shift to multi-user multi-base sta¬ 
tion scenarios, enabling on one hand exploitation 


of multi-user diversity in space, time and fre¬ 
quency and, on the other hand, consideration of 
interference in-between simultaneous transmis¬ 
sions from multiple base stations. 

We will address cross-layer multi-user schedul¬ 
ing, jointly optimizing multi-user resource allo¬ 
cation and per-user link adaptation; this is an 
intricate issue in LTE, due to the non-linear rela¬ 
tionship between the resources assigned to a user 
and its corresponding SC-EDM SINR ([5]); we 
have already addressed this issue for the down¬ 
link in m- Multi-user scheduling, furthermore, 
has to find a favourable trade-off between trans¬ 
mission efficiency and fairness of resource alloca¬ 
tion. We will extend existing downlink sched¬ 
ulers, which enable Pareto-efficient transmission 
with arbitrary fairness, to the uplink specihcs 
and compare to other proposals, e.g., m- 

The realization of massive MIMO in LTE com¬ 
pliant systems is another highly important re¬ 
search topic, since it promises an order of mag¬ 
nitude network efficiency gains through spatial 
multiplexing of users |52l ]53[ [Si]. Yet, many 
issues still need to be better understood and re¬ 
solved to enable efficient massive MIMO trans¬ 
mission in practice. One important step towards 
reasonable performance investigation of massive 
antenna arrays is to employ realistic channel 
models, such as, the 3GPP three-dimensional 
channel model m, which we plan to incorpo¬ 
rate in future releases of our simulator. 

We hnally plan to implement multi-base sta¬ 
tion support in future releases of the Vienna 
LTE-A Uplink Simulator. Even though, for rea¬ 
sons of computational complexity, simulations 
will be confined to comparatively small scenarios 
containing some few base stations, we still expect 
to extract valuable performance indicators for 
coordinated multipoint reception schemes [56], 
accounting for practical constraints, such as, lim- 
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30dB SNR. 


Figure 12: Comparison of channel interpolation techniques using different numbers of reference 
symbols and LS-SAV channel estimation. 


ited back-haul capacity. 

9 Conclusion 

For an LTE-A uplink transmission model we de¬ 
rived SINK expressions, both with and without 
DFT pre-spreading. We specialized these equa¬ 
tions to ZF and MMSE receivers and showed 
that ZF performance is strongly affected by the 
worst subcarrier. Comparing the resulting BER 
we revealed, that SC-EDM performance is gener¬ 
ally inferior to OFDM and that applying MMSE 
equalization is crucial to get closer to OFDM 
performance. 

Based on the system’s SINK we analysed the 
achievable rate. We also introduced a method to 
estimate the SC-FDM rate for AIr > L. Further 
a possible calculation of LTE-A link adaptation 
parameters was proposed to achieve throughout 
close to performance bounds. 

Lastly, we considered methods to gather CSI 


at the receiver. We compared the performance 
of various channel estimation and interpolation 
techniques. By incorporating the channel esti¬ 
mates of the previous subframe, we showed su¬ 
perior performance in terms of channel interpo¬ 
lation. 

Appendix 

General MIMO SC-FDM A SINK ex¬ 
pression 

The signal estimates are described via the input- 
output relationship Equation ([I]). We first slice 
out that part of K which acts on layer I by mul¬ 
tiplying with the selector matrix from left. 
As indicated in ([T]), the signal estimate consists 
of three contributions. 

signal: Xg = {I © K) x 

interference: x\ = {K — I Q K) x 
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noise: = n 

As X and n are zero mean random quantities, 
their power is described by means of the second 
moment. To calculate the second moments we 
take out the diagonal elements of the respective 
covariance matrices of each contribution. 


SINRfC-FDM ^ 


(43) 


O O + / 0 Ejsn®^}) 


-1 


1,1 


Before we derive the different covariance ma¬ 
trices, we recapitulate a required property of 
circulant matrices. A circulant matrix C G 

f^NxN 

is fully described by its first column c, 
as its eigenvectors are the DFT basis-vectors 
and its eigenvalues are the DFT of c = 
(cq, Cl, ... , CAr_i). 

/ Co CAT-l ... Cl \ 


C 


Cl 


Co 


C2 


(44) 


\ CAT-i ... Cl Co / 


= Diag {Dc) D = AD (45) 


The main diagonal elements cq of C are given by 
iV-l N-l 

(46) 


Ejisa;^}: 

The input-output matrix K is of block-circulant 
structure, as illustrated in Figure [3] (c). The 
eigenvalues of the diagonal blocks are given by 


diag (A) = 5*^^^diag {FHes) and the diago¬ 
nal elements of the diagonal block are then 
^I'^S'^^Miag (FiTeff) as asserted by Equation 
(I46]l . thus 


SW (10 K) = ll^5Wdiag(FF/efr) I (47) 


Assuming zero mean,white data with variance cr^ 
the diagonal elements of Mix^xl^} are given by 


E{a:ia:(^}: 

If C is circulant 


C = C-coI 


( ° 

Cl 


C7V-1 

0 


Cl ^ 
C2 


\ CN-l ... Cl 0 / 

(48) 

is circulant as well and the diagonal elements of 
CC are the sum of the magnitude squares of 
c = (0, Cl,..., ctv-i). Using Parseval’s theorem 
we arrive at 


N-l 



J_ s^N—1 
N 2^j=l 


I [Mj,j 


(49) 


J_ ^N—l 
N Ajj=0 


I [A] 


h3\ 



•The inter-layer interference consists of L — 1 C- 
type blocks, where we simply average the mag¬ 
nitude squares of the eigenvalues, i.e., the corre¬ 
sponding block-part of FH^^. The intra-layer 
interference is described via a C block and is 
given in Equation (j49h . Both contributions can 
be compactly written as 

(50) 
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Ejainsf}: 

The noise covariance matrix is circulant as well 
and the detailed derivations can be found in [2B] . 

SISO MMSE SC-FDMA SINK expres¬ 
sion 

For a SISO system and an one-tap equalizer 
the expression FHes is of diagonal shape. [26] 
has shown, that the MMSE equalizer for SC- 
FDM equals the OFDM expression, i.e., F = 
(%/ + Thus, the elements on 

the main diagonal of FHgf[ are simply given by 
and we rewrite ([5|) to (l5H) . 
where we have used the identity 




^sc 

E- 


\Hk\^ 




= 1 - 


1 


NsC 

Nsc 




(51) 


from [H]. 
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